Solitons on the edge of a two-dimensional electron system 



OO 

o 

Q 



C. Wexler and Alan T. Dorsey 
Department of Physics, University of Florida, Gainesville, Florida 32611-8440 
(1 Oct 1998, revised 5 Dec 1998) 

We present a study of the excitations of the edge of a two-dimensional electron droplet in a magnetic 
field in terms of a contour dynamics formalism. We find that, beyond the usual linear approximation, 
the non-linear analysis yields soliton solutions which correspond to uniformly rotating shapes. These 
modes are found from a perturbative treatment of a non-linear eigenvalue problem, and as solutions 
to a modified Korteweg-de Vries equation resulting from a local induction approximation to the 
nonlocal contour dynamics. We discuss applications to the edge modes in the quantum Hall effect. 

PACS numbers: 73.40.Hm, 02.40.Ma, 03.40. Gc, ll.10.Lm 
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Shape deformations are important for our understand- 
ing of such diverse problems as the low-lying excita- 
tions of atomic nuclei [jjj, the hydrodynamics of vortex 
patches , the evolution of atmospheric plasma clouds 
j| , the formation of patterns in magnetic fluids and su- 
perconductors H , and the dynamics of suspended liquid 
droplets to name just a few examples. In many 

of these systems there is a field which is approximately 
piecewise constant (for example, the vorticity for the vor- 
tex patches, or the magnetization for the magnetic flu- 
ids); this greatly simplifies the study of the dynamics, 
since it is often possible to focus attention on the bound- 
aries at which this field is discontinuous, and construct 
the contour dynamics for these boundaries, subject to 
some global constraint (such as area conservation) . 

The edges of a two-dimensional electron system 
(2DES), and in particular the edges of a quantum Hall 
(QH) liquid, present a unique opportunity to study the 
dynamics of shape deformations in a clean and controlled 
environment. The 2DES in the QH state is incompress- 
ible, so that the electron density is approximately piece- 
wise constant, suggesting that a contour dynamics ap- 
proach to studying the droplet excitations is viable. In 
addition, the charged nature of the system facilitates the 
excitation and detection of deformations of the droplet. 

In this paper we will formulate the study of the exci- 
tations of a droplet in a 2DES as a problem in contour 
dynamics. In the usual treatment of the edge excitations 
H, a linearization of the equation of motion is done at 
early stages, thus limiting the applicability to small de- 
formations of the edge of the system from an unperturbed 
state. In this Letter we consider non-linear terms which 
are present in the full contour dynamics treatment. We 
first present perturbative results for non-linear deforma- 
tions of the 2DES shape. We then show that the lo- 
cal induction approximation to the full contour dynam- 
ics generates the modified Korteweg-de Vries (mKdV) H 
equation for the curvature dynamics; the mKdV equation 
also arises in studies of vortex patches |3| and suspended 
liquid droplets pj . The mKdV dynamics conserve an in- 
finite number of quantities, including the area, center of 



mass, and angular momentum of the droplet |p^ , |Tl| , so 
that our local approximation to the nonlocal dynamics 
preserves the important conservation laws. The mKdV 
equation also possesses soliton solutions, including trav- 
eling wave solutions which represent uniformly rotating 
deformed droplets. 

/. Hydrodynamics of a two-dimensional electron system 
in a magnetic field. — Consider a 2DES in a strong mag- 
netic field. The electron configuration can, in most cases, 
be characterized by its density n(r) and its velocity v(r). 
Treated as a classical fluid, the dynamics is determined 
by the Euler and continuity equations, 
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where uj c = eB /m e is the cyclotron frequency and e is the 
dielectric constant of the medium; the first term on the 
right hand side of Eq. ([[]) is the Lorentz force, the second 
is the Coulomb interaction, and E cxt is the electric field 
due to the background positive charge, gates, etc. It is 
noteworthy that this simple hydrodynamic theory of edge 
deformations, and the corresponding canonical quantiza- 
tion of the classical Lagrangian of the fluid, is also able 
to capture the essence of the many-body problem in both 
integer and fractional QH states [0. 

The theory of small deformations of the edge has been 
extensively studied in Refs. Jl^-|lJ|. The main conclu- 
sion is that for strong magnetic fields, when Landau 
level quantization becomes important, the only low en- 
ergy modes are edge modes which propagate in one direc- 
tion along the edge of the 2DES; the bulk modes become 
gaped with a minimum frequency lo c . For a sharp elec- 
tron density profile, this edge mode is the "conventional" 
edge magnetoplasmon mode, with a dispersion relation: 
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where k is the mode wave-number, 7 sa 0.5772 is the Eu- 
ler constant, and a is a short-distance cut-off Q (the 
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largest of the transverse width of the 2DES, the magnetic 
length, or the width of the compressible edge-channel). 
In addition, for the more realistic wide and compress- 
ible edges, new branches of "acoustic" modes exist [ p^[ . 
Recent time-of-flight measurements pl| on 2DES in hct- 
erostructures confirm this simple picture. 

77. Edge modes: dynamics and kinematics. — For what 
follows, let us consider the simplest case: a sharp edge 
between an incompressible charged fluid and "vacuum" 
(Fig. |l|) . We focus on this mode since it has the simplest 
structure and is the most readily observable |l2|,|l5|]. As 
in Ref . fl2| , we neglect inertial terms on the left hand side 
of Eq. ([j])^ and find that the electron velocity is given by 
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where v cxt = — (e/m e uj c ) e z x E oxt is the velocity induced 
by the external field, and A is the area of the droplet. 
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FIG. 1. A charged incompressible liquid in a magnetic 
field. We assume a piecewise constant electron density (n = n 
inside, while n — outside the droplet). The parameteriza- 
tion R(ip,t), the tangent t and normal n unit vectors to the 
boundary V are indicated; s is the arclength and 8 the tangent 
angle. 

Let us now concentrate on the "internal" velocity given 
by the first term in Eq. (||) jTJ]]. For an incompressible 
2DES with a piecewise constant density (l7) , the density 
can be taken outside the integral; then using Stokes' the- 
orem, the area integral can be transformed into a line 
integral over the boundary T~dA of the electron liquid: 



v(r) = 
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Here t(s') = tVr(s') = r s / is the unit tangent vector at 
the arc-length s'; for later use, we also define nit as 
the unit normal vector. The short distance singularity in 
the integrand is cut off at a length scale ro . Equation (|5|) 
forms the basis of our contour dynamics treatment — it 
expresses the velocity of the edge in terms of a nonlocal 
self-interaction of the edge. 

For the sake of comparison, we draw analogy to the 
case of a vortex patch, a two-dimensional, bounded re- 
gion of constant vorticity lo p surrounded by an irrota- 
tional fluid, where the interaction is logarithmic 
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The vorticity and the area of the patch are conserved, 
a consequence of Kelvin's circulation theorem We 
see that in both cases (i) the dynamics is chiral, being 
determined by the tangent vector; (ii) the fluid contained 
within r is incompressible, so that the area is conserved. 
It is this analogy which inspired the present work. 

Having determined the velocity of the electron liquid, 
we now focus on the motion of the 2DES boundary F. 
The velocity of a point on the boundary can be written 
in terms of the normal and tangential components, 



v = U(s) n + W(s) t. 
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The tangential velocity W(s) is largely irrelevant. We 
now ask whether there are modes which propagate along 
the boundary with no change in shape. Previous work 
]l2[ has focused on small perturbations of a straight, in- 
finite edge. Here we consider deformations of a circular 
droplet of incompressible electrons. A uniformly prop- 
agating mode is therefore characterized by a boundary 
that moves like a rotating rigid body, namely the radius 
of the boundary satisfies R(ip , t) = R(ip — ilt) , where (f is 
the azimuthal angle and f2 is the angular frequency. This 
translates into a condition for the normal velocity: 

[/ee n(r).v(r)| rer = f!n(r).(e,xr) . (8) 

We seek boundary shapes R((f) (see Fig. [j]) that ro- 
tate uniformly, satisfying Eq. (||) . Consider the following 
parameterization of the surface: 
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The unit tangent vector is given by t((p) = t{lp) /\t(lp)\ 
wheie T(ip) = dr(ip)/dip = e r R'(ip)+e ip R(ip). Likewise, the 
unit normal vector is given by n((^) = — e z x T(ip)/\r(ip)\. 
The normal velocity on the boundary can be written as 
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It has not been possible to solve exactly the non-linear 
eigenvalue problem for bi and Q [Eqs. j8]-^0|)] exactly. We 
therefore seek a perturbative solution by expanding the 
right-hand side of Eq. (JlTJ) in powers of bi. This allows 
us to to go beyond the linear approximations used in the 
past, and we have succeeded in calculating shape defor- 
mations to 0[bf] and angular frequencies to 0[fef]. Ex- 
panding the non-linear eigenvalue problem to fifth order, 
we find the condition 

^ ( bl + \ E b i-v h P ) = Qi b i + ^2 Rl 'P b i-p h P 

V p ) p 
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where fl = (em e oj c Ro/ne 2 ) fl, and the "matrix elements" 
Q, R, S, T and U are obtained from Eq. ([!(]). The details 
of this calculation will be published elsewhere . The 
zeroth order angular frequency for the eigenmode with 
rotational symmetry C; (rotations by 2tt/1) is 
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where the last sum can be related to the digamma func- 
tion ip(\l\ + l/2). This linear result has been previously 
derived by several authors O]; corrections are 0[bf}. For 
a direct comparison with Eq. (|^), which corresponds to 
the large-l limit, we substitute the asymptotic expansion 
for the sum in Eq. (|l2|); multiplication by Rq yields the 
propagation velocity: 
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which closely corresponds to v g = duio(k)/dk after the 
substitution l^ka [see Eq. (0)]. The dispersion for these 
linear edge excitations has been confirmed experimen- 
tally in both the frequency [|l9| and time Jl^,^] domains. 

Some of the shapes are shown by the dotted curves 
in Fig. ||, and closely resemble the "V states" of vortex 
patches found by Deem and Zabusky (2l]]. For larger 
deformations, the appearance of oscillations indicate that 
higher order terms are needed, since they correspond to 
higher Fourier components. 

777. Local induction approximation. — As we have seen, 
the motion of the edge is determined by the velocity of 
the fluid at the surface. The nonlocal equation for the 
velocity of the boundary, Eq. (^|), can be turned into a 
differential equation for the curvature of the boundary 
if we concentrate on the local contributions. This lo- 
cal induction approximation (LIA) jl8| was explored by 
Goldstein and Petrich in a series of papers dealing with 
the evolution of vortex patches P,pj|. The situation is 
considerably more favorable in this problem, due to the 
more rapid decay of the interactionjl / r for charges vs. 
ln(r) for vortices, see Eqs. (||) and (0)]. 

The LIA is an expansion of the integrand of Eq. (||) 
about s' = s; a long-distance cut-off A is introduced by 
replacing the line integral ^ r {- • -}ds' by J/J^wjf {• ■ -}ds'. 
Using the Frenet-Serret relations r s =t, r ss =t s = — kth, 
where k = S is the local curvature of the boundary, to 
lowest order the normal and tangential velocities are p0[ : 
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Since the rate of change of the area of the droplet is 
At — § dsU(s), the LIA, with L/lia oc k s , automatically 
conserves area. 



The time evolution of a curve in two dimensions is 
given quite generally by the differential equation pjlB2[ 



= -[n 2 + d ss ]U + k s W - k s [kU + W sl ]ds' . (15) 



Introducing the results from Eq. ( |14[) , changing "gauge" 
(by modifying W) and rescaling time, we find that the 
curvature satisfies the mKdV equation M: 



Kt — — K K s T" K sss . 
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The mKdV dynamics are integrable, with an infinite 
number of globally conserved geometric quantities 
the most important of which are the center of mass, area, 
and angular momentum of the droplet. 

The mKdV equation possesses a variety of soliton solu- 
tions, including traveling wave solutions and propagating 
"breather" solitons. Here we will focus on the traveling 
wave solutions of Eq. ( |l6| ) of the form n(s, t) =g{z) with 
z = s — ct, which represent uniformly rotating deformed 
droplets. The ordinary differential equation for (7(2) can 
be integrated twice with the result 
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where a and b are constants of integration (a = 6 = for 
infinite systems). The periodic solutions of this equation, 
expressed in terms of Jacobi elliptic functions, are: 
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where K max ,min are the maximum and minimum cur- 
vatures, and £ is determined by the boundary condi- 
tions. The period of k is given by the elliptic integral 
L K = (8/y/pq) K{\). The first integral of the curvature 
gives the tangent angle 9{s) — J* n{s')ds' . We require 
that 8(1 L K ) = 2ir, so that the resulting curve is closed 
(this fixes £). The integer / determines that the curve 
has Ci symmetry. The curves thus generated can be char- 
acterized by ((,K m ax,^min) or more conveniently by the 
symmetry, the area and the perimeter of the curve. 

The full lines in Fig. ^ show some uniformly rotating 
soliton shapes, calculated from Eq. (|T8|). These are es- 
sentially identical with Goldstein and Pctrich's |ll[ soli- 
ton solutions for the vortex patch problem (Deem and 
Zabusky 's "V states" ^lj), but the more local interac- 
tion in the QH case guarantees a better correspondence 
with the perturbative solutions. Indeed, the curves re- 
sulting from the perturbative method and the LIA are 
quite close, even for considerable deformations of the 
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boundary. For larger deformations the perturbative re- 
sults show artifacts due to the limited number of Fourier 
components. The advantage of the LIA becomes evident 
in this case, since it is an expansion in powers of the 
curvature and not the deformation, and thus allows for 
relatively large long- wavelength deformations. More sig- 
nificantly, the LIA and the resulting integrable dynamics 
allow one to uncover geometrical conservation laws which 
would be hidden in a perturbative calculation [|l0|. This 
advantage comes at a price: the detailed information on 
frequencies is obscured by the introduction of the long 
distance cut-off A and by the gauge transformation of 
the tangential velocity W, while the frequency is easily 
obtained in the perturbative calculation. 
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FIG. 2. Uniformly rotating shapes of a 2DES. Solid lines: 
solutions of the mKdV equation obtained from the local in- 
duction approximation. Dotted lines: solutions obtained us- 
ing the perturbative expansion. The values of the coeffi- 
cient bi, and the ratio of curvatures a = K m in/«inax are: (a) 

62 = 0.073, cr = 0.4; (b) 6 2 = 0.19, (J = 0; (c) fe 2 =0.36, a = -0.2; 
(d) b 3 = 0.027, a = 0.4; (e) b 3 = 0.10, a = -0.2; (f) 

63 = 0.29, a = -0.45; (g) 64 = 0.014, a = 0.4; (h) 

64 = 0.089, cr = -0.4; (i) 64 = 0.24, a = -0.56. 

IV. Conclusions. — A contour dynamics formulation of 
the excitations on the edge of a two-dimensional elec- 
tron system in a magnetic field has allowed us to demon- 
strate the existence, beyond the usual linear regime, of 
shape deformations that propagate uniformly. A local 
approximation to the nonlocal dynamics shows that the 
curvature of the edge of the droplet obeys the modified 
Korteweg-de Vries equation, which has integrable dy- 
namics and soliton solutions. Since these solutions are 
dispersionless, it may be possible to distinguish them 
from linear edge waves in time-of-flight measurements. 
Earlier studies |^] have shown non-linear waves in edge 
channels, but the origin of the non-linearity and exper- 
imental details are different than considered here. On 
the theoretical side, it would be interesting to connect 
our hydrodynamic treatment of these edge solitons with 
field-theoretical treatments of edge excitations O] . 
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